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Electromagnetically induced transparency and absorption (EIT and EIA) are usually demonstrated using 
three-level atomic systems. In contrast to the usual case, we theoretically study the EIT and EIA in an equivalent 
three-level system: a superconducting two-level system (qubit) dressed by a single-mode cavity field. In this 
equivalent system, we find that both the EIT and the EIA can be tuned by controlling the level-spacing of the 
superconducting qubit and hence controlling the dressed system. This tunability is due to the dressed relaxation 
and dephasing rates which vary parametrically with the level-spacing of the original qubit and thus affect the 
transition properties of the dressed qubit and the susceptibility. These dressed relaxation and dephasing rates 
characterize the reaction of the dressed qubit to an incident probe field. Using recent experimental data on su- 
perconducting qubits (charge, phase, and flux qubits) to demonstrate our approach, we show the possibility of 
experimentally realizing this proposal. 



PACS numbers: 85.25-j, 42.50.Gy, 42.65.An 

I. INTRODUCTION 

A. Electromagnetically induced transparency in optics and 
superconducting circuits 

Electromagnetically induced transparency (EIT) 01 0] 
manifests spectroscopically the quantized three-level structure 
of an atomic medium through its interactions with two semi- 
classical fields. It has been widely explored in various con- 
texts (e.g., [3-6]) since its inception. For example, the EIT 
effect has been studied in the context of a two-level atom 
instead of the usual three-level system. In Ref. 0], the en- 
ergy levels of the atom are split by a driving optical field into 
doublets, equivalent to ac-Stark shifts, and the transparency is 
realized on the final four-level system. 

The development of superconducting quantum circuits 
(SQC) in recent years have heavily employed concepts from 
quantum optics, and SQC have become a testbed of quan- 
tum optical phenomena, including EIT (e.g., [7-9]). In one 
case 13], a flux qubit system has even taken advantage of the 
EIT effect as a means to probe decoherence by circling the 
flux qubit with a readout-SQUID in the SQC. Both the optical 
version and the SQC version of the EIT effect are illustrated 
on the left and right panels, respectively, of Fig.Q] to compare 
their similarity and differences. This also serves as a prelude 
to our discussion of the tunability of electromagnetically in- 
duced transparency and absorption. 

For the optical version, on the left of Fig. [JJ two different 
kinds of classical electromagnetic waves O p (t) = f2 p e~ ? " pt 
and fi c (t) = il c cxp _l " c *, are shown, for a typically weak 
probe field f2 p (i) and a strong control field f2 c (t) (where 
f2 c <C f^ P ). These are shown incident on an optical medium in 
Fig. 02a) and (b), respectively. The frequencies, uj p and u c , of 
these two fields are resonant with some energy level spacings 
of the optical medium, which are typically A-type three-level 
atoms, and each of these two fields alone will be absorbed 



and cannot travel through the medium. Figure |TJg) shows a 
canonical level diagram for such three-level atoms, where the 
probe field is resonant with the level spacing between |1) and 
1 3) and the control field with that between |2) and |3). 

When the probe and control fields are simultaneously in- 
cident on the medium, as shown in Fig. [TJc), the resonance 
between the control field and the medium will render the 
medium detuned from the probe field and hence let the probe 
field travel through without being absorbed. The level |3) 
is driven out of its original position from |2) and make the 
1 1 ) -to- 1 3) spacing detuned from the probe field frequency in 
Fig. |TJg). In other words, an absorbing medium becomes 
transparent to an incident probe field when a control field is 
simultaneously applied. More precisely, the transparency of 
the probe field can be considered as an effect of its quantum 
interference with the control field. The strong coupling of 
the medium with the control field perturbs the original level 
spacings and provides two excitation pathways of equal prob- 
ability but opposite signs to the probe field, indicated by the 
Autler-Townes levels |3+) and |3— ) in Fig. [Tig)- The result- 
ing signal that exits from the medium is thus a destructive 
superposition of two versions of the same signal, hence a de- 
structive interference and a zero-absorption of the linear sus- 
ceptibility (see, e.g., 0] for a comprehensive review). Note 
that in the atypical case where both the probe and the control 
are strongly coupled to the medium, the interference pattern 
is severely altered: in one case, simultaneous transparency for 
both fields is achieved lTT"oh whereas in another, an enhanced 
absorption of the probe can occur ifTTll . 

The SQC version of the optical EIT is illustrated on the 
right side of Fig. [T| where the optical medium is replaced by 
a Josephson-junction multi-level system. In a similar man- 
ner, when the probe and the control signals np n (t) and OJ, n (i), 
which are either current or voltage signals in this case, are fed 
separately into the SQC [shown in Figs. Old) and (e)], the am- 
plitudes of their outputs £7° ut (i) and fi° ut (t) will be smaller 
than those of their inputs. This shows that the SQC sepa- 
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FIG. 1 : (Color online) Schematic diagram illustrating electromagnet- 
ically induced transparency, or EIT. Panel (a) schematically shows a 
medium that strongly absorbs a probe light beam of amplitude Q p ; 
(b) shows the same absorption, but for a control light beam of am- 
plitude fi c . However, when both beams are applied simultaneously, 
as in (c), then the previously absorbing medium becomes transparent 
for the probe beam Q, p , i.e. transparency is electromagnetically in- 
duced by the control field fi c . The superconducting circuit analog of 
this phenomenon is schematically shown in (d), (e), and (f), where a 
"box" schematically represents the circuit. In (d) a classical voltage 
or current signal (in the microwave range) of amplitude is fed 
into the input of the circuit. The signal is "absorbed" in the circuit 
and the amplitude Q p ut of the output signal is less than that of the in- 
put, (e) shows the same effect with another microwave control signal 
of amplitude fi c . However, when both microwave signals are applied 
simultaneously to the circuit, as in (f), then the previously "absorbing 
medium", now a circuit, becomes transparent for the first signal Q,' p . 
(g) shows a canonical level diagram for a three-level system coupled 
to a weak probe field and a strong control field. 



rately absorbs the input signals due to resonances between the 
input electrical signals and the energy level spacings of the 
Josephson-junction device. Nonetheless, a simultaneous ap- 
plication of the two input signals would let the probe signal 
conduct through the circuit without losing energy, while the 
control signal is mostly absorbed, as shown in Fig.[TJf)- 



B. Electromagnetically induced absorption 

Two closely-related, but far less studied, optical phenom- 
ena are electromagnetically induced absorption (EIA) lfl3ll 
and switchable dispersion lfl411 . where the hyperfine structure 
of the ground state of an atom is used. The quasi- or near- 
degenerate levels originate from the same ground state hyper- 
fine level and have a very small splitting; when the configu- 
rations of their total angular momentum and that of their ex- 
cited state have even parity, a A-type three-level system with 
closed cyclic transitions is formed. Contrary to the odd-parity 
case, where the three-level system becomes A-type and can 
exhibit EIT, the even-parity configuration makes the multi- 
level medium absorptive to the probe field even when it is 
resonant with the control field. These relations between opti- 
cal properties and the parities of SQC were recently predicted 
theoretically lfl5ll and verified experimentally iflql based on 
selection rules and symmetry-breaking. 



C. Tunable transparency and absorption 

The circuit designs in Ref. 10, 01 are based on the multi- 
level energy structure of Josephson junction devices. These 
designs require advanced measurement techniques where the 
third and higher energy levels are often far separated from the 
bottom two. We therefore consider in this article an alterna- 
tive way to construct a multi-level energy structure in SQC 
by "mixing" a two-level system with a resonant field, form- 
ing a dressed multi-level structure that is tunable and does not 
entirely rely on the device characteristics of the junctions. 

Considering the recent progress in studies on superconduct- 
ing qubits (e.g. lfl7l - l2lll ) and dressed superconducting qubits 
by a single-mode cavity field I22H25I1 . we study here a dressed 
three-level qubit equivalent to those in three-level Josephson 
junction devices, in order to theoretically realize the effect of 
EIT on SQC's. Also motivated by the studies of EIT and EIA 
in atomic systems (e.g., filial [l3l [Till ), we show how the EIT 
and the EIA phenomena coexist and transmute into each other 
on the same dressed SQC. 

We will select three energy levels among the multiple lev- 
els of the dressed superconducting qubit. The tunable level 
spacing of the superconducting qubit then not only affects the 
level splitting of the dressed states, but also affects the relax- 
ation and dephasing rates of the dressed system. These tun- 
able relaxation and dephasing rates will determine the sys- 
tem's specific dynamics when coupled to a classical signal 
field, effectively making it have either a A-type (closed) or 
a A-type (open) transition pattern. Therefore, if two classical 
electromagnetic fields, the probe field and the control field, 
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are fed concurrently into the circuit, the complex susceptibil- 
ity of the dressed qubit gives an absorption spectrum that ei- 
ther dips or peaks at the zero probe field detuning. The choice 
of dip or peak of the spectrum is analytically determined by 
a biquadratic equation which is dependent on the qubit level 
spacing and the environment temperature through the dress- 
ing process. By determining the number of real roots given 
by this equation, we can distinguish two effective regimes of 
operations: EIT and EIA. 

The dip-type and the peak-type spectra corresponds to the 
EIT and the EIA effect, respectively, depending on the mag- 
nitude of the qubit level spacing with respect to the eigen- 
frequency of the resonator quantum field. Compared to the 
atomic case, the lower two dressed states from the super- 
conducting qubit hence act like those hyperfine levels from 
the atomic ground state, giving either a closed transition for 
the EIA regime and an open transition for the EIT regime. 
Note also that these tunable regimes of operations are closely 
related to the tunable luminosity suggested by Agarwal et 
al. H2I1 . where the group velocity of light is controlled by a 
"knob" signal field that couples the metastable states. The tun- 
able metastable-state coupling given effectively by our dress- 
ing process is thus comparable to that given directly through 
the knob field. 

Our analysis will focus on superconducting quantum cir- 
cuits that have a strong coupling between the qubit and the 
resonator l23l l24ll . realized by a coplanar waveguide (CPW) 
transmission line. These circuits include combinations of the 
CPW resonator with either a charge qubit I25U27I1 . a phase 
qubit ll28M30h , or a flux qubit (HI. 

In Sec. m we first describe a general theoretical model and 
derive the energy spectrum of the dressed multi-level system. 
The dressing process is described in Sec. [HI] The first-order 
susceptibility and the dressed relaxation and dephasing rates 
among the multiple levels of the qubit-resonator combination 
are calculated in Sec. [IV] using a density matrix formulation. 
The determination of the switching between the transparency 
and the absorption as well as the discussion of the correspond- 
ing transition patterns are presented in Sec. [VI] In Sec. lVIIl we 
then consider experimentally accessible parameters for differ- 
ent types of qubits to demonstrate and numerically analyze our 
theoretical results. Conclusions are summarized in Sec. I Villi 



II. UNDRESSED QUBIT 

Our discussion of the superconducting qubit system is in- 
dependent of the specific type of qubits (phase qubit, charge 
qubit, or flux qubit) employed. Because the two-level struc- 
tures of these qubits are commonly described by a general 
Hamiltonian with a a z term and a a x term ltr7H20ll . To sim- 
plify our discussion, here we consider the qubit in the diagonal 
basis whose eigenfrequency u> q is the root-mean-square of the 
coefficients of the a z and a x terms. The CPW resonator that 
couples to the qubit, akin to a cavity for photons, is described 
by a pair of annihilation and creation operators a and with a 
resonant frequency denoted by loq. The dipole-field coupling 
between the qubit, acting as the dipole, and the energy quan- 
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FIG. 2: Schematic diagram of the "dressing" process: the Fock num- 
ber states |n) of a coplanar- waveguide resonator (shown in the first 
column) as well as the ground state \g) and excited state e) of a 
qubit provide the "undressed" tensor-product states jn, g) and \n, e) 
(shown on the middle column). These tensor states are modified by 
their mutual interaction, producing the renormalized or "dressed" 
states |/i n ) and \i/ n ) (shown on the right). Mathematically, these 
states correspond to the eigenvectors of the circuit Hamiltonian in a 
non-diagonal basis (undressed states) and the diagonal basis (dressed 
states), respectively. 



rum within the CPW resonator, acting as the cavity field, is 
along the z-direction of the qubit in the non-diagonal basis. 
In the diagonal basis of the qubit, the Hamiltonian between 
the qubit and the transmission line resonator is given by the 
Jaynes-Cummings model 



(1) 



with a rotating-wave approximation and H = 1, where 77 de- 
notes the coupling constant between the cavity field and the 
qubit. Here, for convenience, we still use a z to denote the 
qubit operator in the diagonal basis. 

Two signals, the probe signal and the control signal, with 
traveling frequencies denoted by uj p and uj c and Rabi frequen- 
cies denoted by fl p and il c , respectively, are fed into the cir- 
cuit. These two signals are treated as classical electromagnetic 
fields and their interaction Hamiltonian with the qubit can be 
written as 

#ext = Q p e iu;p V_ + O c e iWc V_ + h.c. (2) 

The total system Hamiltonian is then given by 

H = H cir + ifext- (3) 

The energy states of the CPW resonator and the qubit before 
dressing are shown as horizontal lines on the left part of the 
schematic diagram in Fig. [2] Note that the energy levels \n,g) 
and |n, e) have equal spacings for all n: for fixed probe field 
and control field, the undressed qubit with a particular level 
spacing might not resonate with them. Having been dressed 
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by the CPW resonator, the qubit will exhibit a spectrum with 
numerous energy levels spaced in a tunable non-uniform pat- 
tern (the dressing process will be discussed in next section), as 
shown on the right part in Fig. [2 providing more possibilities 
for matching levels between the dressed qubit and the probe 
and control signals. 

III. DRESSED QUBIT 

The dressed qubit Hamiltonian can be derived by rewrit- 
ing the Jaynes-Cummings model, which usually describes 
the atom-photon coupling for the circuit Hamiltonian H C - 1T in 
Eq. (0. Note that the set {|n, e) , \n + l,g)} spans an invari- 
ant subspace V n of H^, where n denotes the number of en- 
ergy quanta in the CPW resonator; while e and g denote, re- 
spectively, the excited and the ground state of the supercon- 
ducting qubit. Therefore, the corresponding Hilbert space, 
that the Hamiltonian H ca in Eq. acts on, can be written 
as the direct sum 

V = {\0,g)}®V n , 

71 = 

where the ground state (of the combination of the qubit and 
the CPW resonator) does not belong to any invariant subspace. 
The basis of each subspace V n can be transformed by rotating 
an angle 

°n = U^(^±L) (4) 

2 \uj q -uj /2J 

such that the circuit Hamiltonian H cn is diagonalized in the 
invariant subspace V n with eigenvalues 



En = U + -Jwo±y(w,- ^) +rf(n + l), (5) 

where the second term constitutes the Rabi splitting for each 
energy level n. Written in the transformed basis 

\n n ) = cos#„ |n, e) - sin# n |n + l,g) , (6) 
\u n ) = sinOn \n,e) + cos#„ \n+ l,g) (7) 

and neglecting the ground state \0,g), the circuit Hamilto- 
nian can be expressed in its diagonal "dressed" form: 

tfcir = J2 ^ E n I/O (Mn| + K K) (v n \] , (8) 
n 

where E^ (E%), associated with the basis vector \/j, n ) (\v n )), 
corresponds to the plus (minus) sign of the eigenvalue of 
Eq. at the n-th Rabi splitting. The corresponding dressed 
eigenvectors are diagrammatically illustrated as the lower (up- 
per) lines on the right side of Fig . |2] 

Note from Eq. that the parameter oj q provides a means to 
tune the level spacing of the superconducting qubits, through 
externally controlling gate voltages, or magnetic flux, or cur- 
rent lfl7l - l20ll . Thus the dressed qubit can exhibit different re- 
sponses to the probe field signal. 



IV. COMPLEX SUSCEPTIBILITY 
A. Three-level system 

Originally, the undressed qubit has two levels. When this 
two-level system is coupled to a driving CPW resonator, the 
interaction "dresses" the system to have an infinite number of 
states, instead of two. These states are shown in Eqs. ©-(0. 
The multi-level structure of the dressed system gives vast se- 
lections of three-level structures on which the EIT or the EIA 
effect can be demonstrated. Before we select three specific 
levels for our purpose, we shall rewrite the external part H cxt 
of the total Hamiltonian, which we have not discussed so far, 
in the transformed or dressed basis. 

The transformed basis spans the product space of the qubit 
space and the resonator space. Considering this, we write 
the flip-up operator <r + = I <E) |e) (g\ = J2 n \ n ' e ) ( n :Sl as 
the tensor product of two space bases in H cxt . Taking the 
inner products of these basis vectors and those of the trans- 
formed basis, we find that, except for the first off-diagonal 
elements (i.e., (fi n+1 \a + \fi n ), (v n+ x\a + \v n ), {^ n+ x\a + \u n ), 
and (v n +i\a+\fi n )), all the entries (including the diagonal 
ones) of the operator cr + in the new matrix representation are 
zero. The #„-dependent non-zero matrix elements are all real 
and, using Eqs. (0-0, the new representation reads 




- cos6> n+ i sin0 n \fj, n +i) {Hn\+sin8 n+1 cosft 



n 



x \v n +i) {v n \ - sin6», i+ i sin6»„ \v n+1 ) (/x n | 

+ cos 0„+i cos 0„ |Mn+i} (v n \\- (9) 

The derivation above applies equally well to the adjoint c_. 

From Eq. (O, we see that the excitations of the qubit engage 
the nearest set of neighboring levels of the CPW resonator in 
a way that the transition coefficients depend on the device pa- 
rameters of the superconducting qubit. Therefore, from the 
point of view of these dressed qubit levels, the process of en- 
ergy pumping into a resonator through a mediating qubit j29tl 
is a laddering of consecutive level -jumps to the next (n + 1) 
dressed level. 

Substituting Eq. (0 and its adjoint into Eq. (0, and also 
using Eq. (0, one reaches a total Hamiltonian expressed com- 
pletely in the dressed basis. Selecting the three levels with 
lowest eigenenergies according to Eq. 0, i.e. 

{|a'o> , H) , |mi>} , 

a Hamiltonian resembling that of a three-level atom with two 
associated dipole-field interactions is obtained 

H A = E$\iM,)<jio\+Eo N>M+^i |A*i><A*i| 
-f2 p e~ iWpt cos 6»i sin^oljUi) (Mo| 
+fi c e~^ ct cos6'icos6'o \m) (u \ + h.c. (10) 

The three selected levels, whose transitions are coupled to 
the classical probe field of frequency lu p and the control field 
of frequency lu c , are shown on the right side of the level dia- 
gram in Fig. Note that when selecting the levels, the state 
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|0, g), which did not participate in the basis transformation, is 
ignored. The selected states (|/xo), \v$), and correspond 
to the dressed ground state, metastable state, and excited state, 
respectively. Each of the classical fields, ui p and uj c , can drive 
transitions between: (i) the ground, metastable, and excited 
states, as well as (ii) those between levels of higher energies. 
We choose the frequencies of the classical fields such that the 
probe field lu p (the control field uj c ) is near-resonant with the 
transition (/io| (|mi) (^ol) f° r our discussion of the elec- 
tromagnetically induced transparency and absorption effect. 
The two frequencies cj p and lj c can be considered to be far- 
detuned from each other and from other transitions (including 
transitions to higher-energy levels). Therefore, all other types 
of interactions can be neglected in the Hamiltonian ilOi . 

B. Qubit-resonator interaction 

The coupling coefficients of the three-level system, given in 
Eq. ( [Tol l, not only depend on the Rabi frequencies il p and f2 c , 
but also on the rotation angles 9q and 9±. That is, the diago- 
nalizing transformation partially determines the magnitude of 
the probe and control field couplings. 

From Eq. (|4j, we observe that the diagonalizing angles at 
two different CPW-resonator levels obey the general relation 

tan 29 n s/rT+T 
tan 29 m yjm + 1 

If the levels are such that m < n, then 6„ > 9 m in the first 
quadrant. Applying this relation to m = and n = 1, we find 

9 1 = ^tan -1 (V2tan26>o) , 

which is not an everywhere-differentiable function with re- 
spect to #o- There exist non-differentiable points at (Itt/2 + 
7r/4), for integer /. Figure [3] plots the rotation-angle- 
dependent factors (cos 0i sin 0o ) and (cos #i cos $o) m the 
coupling coefficients of Eq. 1101 . as functions of 9q, over one 
period, within which four non-smooth turning points can be 
spotted for each function. 

From the first quadrant (0 < #o < 7r/2) of Fig. [3] we note 
that the coupling amplitude between the dressed qubit and the 
probe field is monotonically increasing while that of the con- 
trol field is monotonically decreasing. The two coupling am- 
plitudes coincide at 9q = ir/4. Therefore, the rotation angle 
can be tuned in such a way that the coupling strengths between 
the dressed qubit and the external signals vary considerably, 
from 6q = 0, where the coupling to the control field is maxi- 
mized, to the opposite limit 9q = tt/2, where the coupling to 
the probe field is maximized. 

Coherent trapping, usually discussed in the context of quan- 
tum optics, also occurs here in this superconducting dressed 
two-level system. If the initial state of the dressed qubit is pre- 
pared with equal populations in the two lower states with no 
phase difference, i.e., |V>(0)) = {\po) + l^o)) /v2, then the 
interactions with the external signals, as described in Eq. ( fTOb . 
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FIG. 3: Coupling amplitudes of the dressed qubit to the probe field 
(denoted by the dashed curve) and to the control field (denoted by 
the solid curve). 



will drive the population at the excited level to be dependent 
only on the coupling Rabi frequencies Sl p and tt c . In the 
dressed system, when the qubit level spacing reaches the exact 
values 

w = V 

q 2 tan ^tan-^fic/Op)] ' 

the population at the upper lever will remain zero, while those 
of the two lower levels stay the same over all time t. The 
phenomenon of "trapping" is then achieved in the sense that 
the dressed system will remain in such a configuration with no 
population in the upper excited level, even though the classical 
signals are continually pumped into the system. 



C. Demonstrating EIT via the complex susceptibility 

Using the standard density matrix formalism ll32tl . in this 
subsection, we demonstrate the EIT effect through the deriva- 
tion of the first-order susceptibility of the dressed qubit. The 
dressed qubit acts as a signal-absorbing medium driven by the 
two external signals lu p and us c , as schematically shown in 
Fig. [2] similar to the three-level atom with the optical fields 
given in Fig. Htg)- We will use the matrix element notation 
Pa/3 = \cn) (f3\, where a and j3 can be one of the symbols /i, v 
for the lower levels |yUo) , l^o) an d 1 f° r the excited level |^i); 
p a p denotes level populations when a = j3, or transition am- 
plitudes otherwise. We shall also use the shorthands for the 
coupling coefficients: 

C p (i) = -Q p e- i ^ t cos9 1 ain9 , 
Co ft) = ^ce^ 'cosflicosflo. 

The matrix elements thence evolves according to 
Schroedinger's equation with respect to the Hamiltonian 
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in Eq. ( fTUb . analogous to the optical Bloch equations, 

Pfifi = -TuPw - KvPiti + KpPfMi, (Ha) 

Pvv = -TvPvv - KcPlv + K*Pvl, (Hb) 

pn = — ripn + iCpPip + Kcpiu - i(pPtii - Klpui, 

(lie) 

Pui = - [i( E o - E i) + 7/*i] Pui ~ K P (Pu ~ Pw) + Kcppv, 

(lid) 

pul = - [i(E% - E%) + Jul] Pul - Ko{P\l ~ Pvv) + KpP^vi 

(He) 

Pnv = - - E%) + 7^] pp V - KpPti + 

(llf) 

In the equations above, we have added phenomenologically 
the relaxation rate (T^, T{) for the level |/io) (l^o); \pi)) 
as well as the dephasing rate 7^1 (7^1, 7^) for the transi- 
tion \/j, ) (fii\ (\i/q) |/xo) (fn|)- Note that the subscripts 
/i, z/, and 1 here refer to the levels \/J,o), \vq} and | /Ui) , re- 
spectively, to simplify the notation. This system of equations 
is homogeneous, which gives a zero steady-state solution, i.e. 
the system's thermal equilibrium state, when the coefficient 
matrix is nondegenerate. Therefore, we can assign, without 
loss of generality, the population under consideration entirely 
to the ground state, i.e. p^J] = 1 and p^v = Pi± = 0, and 
the polarization (p( )) (dipole moment) amongst the three 
energy levels at the zeroth-order expansion to zero value, i.e. 

Pfj.1 — Pvl — Ppv — U. 

Since we are concerned with the dispersion and absorption 
spectrum of the dressed qubit, only the first-order perturba- 
tive expansion of the density matrix elements in Eqs. dl ldb - 
( 11 lfl ) are needed. Substituting the steady-state solution above 
into these two equations and removing the time dependences 
of the coefficients in the rotating frame of reference pty — !> 

p^l exp{-iw p i}, p$ -> p$ cxp {-i (w p + Eft - E?)}, 
we have 



=-{iA + 7m1 )p^l + in c cos 9, cos fl o p00 
— iQp cos 9\ sin 8q , 



(12) 



p$ = - (»A + 7/^)^2 + i«c cos 0! cos fl pW , (13) 



where 



(E?-Eg)-, 



is the detuning between the probe field and the level spacing 
of the ground state and the excited state. The control field 
frequency lj c is assumed to match the level spacing between 
the metastable state \i/q) and the excited state \p,i): w c = E± — 

The strength of the probe field is related to the polarization 
of the medium through the relation 



coX 



Eq. ( fT3b in steady state, we can find the first-order susceptibil- 
ity 

X W =x'+ix" 
decomposed in a real part 

X = ZA A 2 - 7^7^ - (fi c cos 6 X cos 9 ) 2 

+ 7^(^1+7^)] (14) 

and an imaginary part 

X " = Z [ 7m1 A 2 + 7^7^ + 7/1 „(fi c cos 9 1 cos O ) 2 ] (15) 
where the common factor is 

l<V| 



CO 



cos 6»i sin 6»o < A ( 7/lX + 7^)' 



where e denotes the vacuum permittivity and cLi the qubit 
dipole moment. Combining this relation with Eq. ( fT2l ) and 



+ [A 2 - 7^17^ - (fi c cos 6»i cos6> ) : 



Note that the dependences of Eqs. (fT4l i-(fT5]l on the detuning A 
are similar to those of the susceptibilities of regular three-level 
atomic systems. The difference here is the extra dependence 
on the qubit level-spacing ui q through the transformation an- 
gles 9\ and 8q. First, the Rabi frequency (f2 c cos^i cos#o)> 
which describes the coupling strength of the control signal to 
the dressed qubit, has an explicit dependence on the angles 
0\ and 60 ■ Second, the relaxation rates 7^1 and 7^ have im- 
plicit dependences on these angles. These cosinusoidal depen- 
dences will vary the magnitude of every term and will hence 
vary the functional behaviors of the two parts x' an d x" of the 
susceptibility. In the following sections, we will see that these 
dependences are crucial to the tunability of the dressed qubit 
to become either transparent or absorptive. 



V. DECAY RATES BETWEEN THE DRESSED LEVELS 
A. Overview of recent theories 

In this section, we study how the relaxation rates 7^1 and 
7ai „ in Eqs. (IT4l)-(IT5b and other dephasing rates are implicitly 
dependent on the qubit level spacing Lu q and other system pa- 
rameters. In particular, we will determine how these decay 
rates are affected by the dressing process, i.e., after the qubit 
levels are rotated into the new dressed basis. 

The decay of an undressed qubit induced by an environ- 
mental coupling has been extensively studied in Refs. l34ll 
and ll35ll through a first-order perturbation model based on 
the Bloch and Redfield's analysis l36ll on what was called a 
relaxation or distribution matrix. The relaxation process in- 
duced by the environment has many implications on multi- 
level systems, including the lasing action of superconducting 
circuits 13711 . 

The recent article by Wilson et al. HHl is the first study that 
took a step further to deal with the dressing effect on relax- 
ation and dephasing, in which the resonator field couples to 
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the z-direction of the qubit and diagonalizing in a displaced 
basis results in a splitting of levels that has a Bessel-function 
dependence on the quantum number n. The approach here 
with the rotating wave approximation, in contrast, has a linear 
dependence of the splitting on n (Cf. Eq. ([5])). This number 
n determines both the thermal distribution of energy quanta 
in the CPW resonator and the rotation angles of the levels in 
each invariant subspace, for which we expect to see a two- 
fold dependence of the relaxation rates on the level n. More- 
over, the environmental temperature T determines the CPW- 
resonator's thermal distribution, as well as the magnitude of 
the fluctuations induced by the thermal coupling of the non- 
dressed qubit. Therefore, the dressed qubit has a two-fold de- 
pendence on the temperature T. 

We will study each aspect of this dependence separately 
in the following subsections. The direct fluctuations intro- 
duced by the thermal coupling before the dressing process 
will be discussed first. The indirect dependence through the 
resonator's thermal distribution after the dressing process will 
follow, where the implicit dependence on the qubit level spac- 
ing uig is also introduced. 



B. Relaxation and dephasing before dressing 

The two steps of computing the dressed decay rates can 
be separated by first assuming that the system has reached 
a thermal equilibrium and the two-level qubit provides the 
only source of noise (i.e., the resonator has no thermal fluc- 
tuations). Following the standard methodology (e.g. in 
Ref. 113411 ). we describe the thermal environment by a sin- 
gle quantum variable X, which interacts with the qubit along 
all three directions through the Hamiltonian H] = (c x a x + 
c.yOy + c z a z )X. Among the interaction coefficients, those of 
the a x - and a y -directions affect the population decay between 
two levels, i.e. they contribute a total relaxation rate 

r 1 =r l + r-f 

which consists of a "down" relaxation rate and an "up" 
relaxation rate ry. 

n = \cx + ic y \ 2 S x {^)/^ 
r\ = \cx + ic y \ 2 S x {-u)/A. 

The Sx shown in the expressions above indicate the power 
spectrum of the heat bath and its value has an exponential de- 
pendence on the environment temperature T, 



1 

2^ 



where R denotes a nominal resistance. 

The a z -direction of the interaction coefficients contributes 
a part, r v , of the total dephasing rate, which is known as the 
pure dephasing part 

r v = \c z \ 2 S x {Q)/2. 



The interactions along a x - and a y -directions contribute the 
other part, 7*1, of the dephasing, which make the total dephas- 
ing rate of the qubit sum up to 

Tl , 
T2 = -+r v . 



C. Relaxation and dephasing after dressing 

Considering the level mixings due to the unitary transfor- 
mations in Eqs. © and (0, the fluctuations of the dressed 
levels become the superpositions of the fluctuations originated 
from the states of the undressed qubit. 

The density matrix element for the lowest level of the 
dressed qubit, expanded in the undressed basis, reads 

| Mo ) ( Mo | = cos 2 6 |0, e) (0, e\ + sin 2 9 \l,g) (l,g\ 

— cos 9 sin# [|0, e) (1, g\ + h.c] . 

We see that this matrix element has the diagonal parts and the 
off-diagonal parts in the undressed basis. To simplify the for- 
mulation, we assume that the dressed qubit has reached a ther- 
mal equilibrium and ignore the system relaxation due to the 
energy exchange between the CPW resonator and the qubit, 
that is, the off-diagonal parts. Then by tracing out the sub- 
space of the resonator part by assuming its energy quanta has 
a Boltzmann distribution, we can obtain the reduced density 
matrix element 

4, = (1 - e-^°) [cos 2 O |e) (e| + e"^° sin 2 9 \g) (g\] , 

where f3 = l/ksT denotes the inverse temperature. Con- 
sidering that the two levels of the undressed qubit contributes 
equally to the up and down relaxations, we arrive at the re- 
laxation rate T^ of the lowest energy level \/j,q) of the dressed 
qubit 



|c a + ^| 2 SxM (1 _ e -^ 0) [cos 2 6o + e -^ sin2 0o ] 

(16a) 

Using similar steps, we can derive the relaxation rates (r„ and 
T{) of the other energy levels (|z/ ) and \ni)) as well as their 
dephasing rates (7^, >y vl , and 7^): 



-p _ \ c x + * c j/| ' 

4 [ 



- e-^°) [sin 2 6 + e-^° cos 2 6 ] 



(16b) 



x [cos 2 6»i + e- f)ua sin 2 61] , (16c) 

7/xi = ~ |C ^ |2 ^ (0) e-**> (1 - e-fr"°) sin 9 cos U (16d) 

7,1 = |c ' |2 f (0) e-^(l - e-^° ) cos 9 cos 9 1 , ( 16e) 

7„ = lCx + l °f Sx{uj) (1 - er^ f cos^o sin^o- (16f) 
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We note from Eq. ( 116fl l that the dephasing rate 7 M „ between 
the split states \p, ) and |i/ ), which share the same number of 
energy quanta, originates from the relaxation rates of the un- 
dressed qubit levels |e) and \g). Therefore, in transforming 
the basis for diagonalization, the roles of relaxation and de- 
phasing have exchanged. 

The total relaxation rate among the three dressed levels be- 
comes now 



r 



which has a two-fold dependence on the environment temper- 
ature T: through the noise spectrum Sx(u) and through the 
resonator distribution cxp(— /3wq). The total relaxation rate 
r is also tunable by the qubit level spacing u> q through the 
transformation angle 6>o- 



VI. TUNING THE ELECTROMAGNETICALLY INDUCED 
TRANSPARENCY AND ABSORPTION 

A. Local extrema of the susceptibility 

The qubit, together with the CPW resonator, acts as a non- 
linear medium of electron propagation in the quantum circuit. 
This nonlinear medium gives the rotation-angle dependent 
and relaxation-rate dependent responses of both the disper- 
sion and the absorption to the incident probe signal S7 p e~ 4Wpt . 
These two responses are quantified by the real and the imagi- 
nary parts of the susceptibility, respectively, in Eqs. <fT4T> - (fl~5T> - 
Since the relaxation rates 7^1, 7^ themselves depend on the 
rotation angles 6>o, #1 (Cf. Eqs. ( 116db -( [T6fl i). which in turn de- 
pend on the qubit level spacing uj q (Cf. Eq. dU), the spacing 
u! q tunes the dispersion and absorption spectra by controlling 
the Josephson coupling energy. 

The absorption spectrum as a function of the detuning 
A between the probe signal and the dressed level spacing 
(-Eq — Ef ) is particularly interesting with its number of max- 
ima depending on the spacing uj q . The derivative of x" with 
respect to A is the product of A and a quartic expression of 
A. Therefore, the absorption spectrum always takes an ex- 
tremal value at the zero root Ao = 0. The quartic expression 
is actually biquadratic, whose two roots are given by 



A+ =±- 



-7/» [Cc(0) +7^17^ 



1/2 



+ (7mi + 7^)Cc(0)^Cc 2 (0) + 7mi7^ , (17) 



where the other two roots that are associated with the negative 
sign of the second term are omitted since the detuning A can 
only admit real values. The two admissible roots coincide and 
meet the zero root, i.e. A± = Aq = 0, when 



That is, at near resonance tu q w wo/2, when the qubit level 
spacing reaches the critical values 



(V2 + l)F 2 + 2\/2 



F\f (a/2 - l) 2 ^ 2 + 16 + 8V2 



(19) 



where 



F(n c ,T) 



u>o\-3 



\c x +icy\ 2 \c z \*Sx{0)S x {u>) 



is a control field amplitude-dependent and temperature- 
dependent factor. We can also check that the condition, 
(Q c cos 9\ cos $o) 2 7ui7«&' 5* 0, has already guaranteed the 
right hand side of Eq. dl7l i to be greater than zero (See Ap- 
pendix A). 

Consequently, inside the critical range Ac- < w g < Ao,+, 
the only admissible root of dx" /dA occurs at the zero point, 
and from the second-order derivative, it can be seen that this 
root corresponds to a local maximum. In the opposite case, 
when the qubit level spacing u q is tuned outside the critical 
range, two new extrema arise in the absorption spectrum, giv- 
ing a total of three turning points in the absorption curve. Be- 
cause x" is an even function of A, the original peak point at 
zero detuning splits symmetrically about the origin, creating 
two symmetric peaks whose distance (A + — A_) is extended 
when oj q is tuned away from the critical values Ac,± at either 
end; whereas the original peak 



x 



1 + fi 



2 cos 81 cos 6*o 



60 



(Q c cos#i cos6> ) 2 +7^i7a^ = °- 



(18) 



itself, where C denotes some constant, starts to dip from a 
maximum value to a local minimum. This local minimum 
tends to zero when the spacing u q is tuned away from its res- 
onant value wo/2 and the amplitude £l c of the control field is 
increased. 

This "peaking-to-dipping" transition (or the increase of the 
number of extrema) indicates the switching of the dressed 
qubit from being transparent to being absorptive to the probe 
field. The magnitude of the qubit level spacing oj q with respect 
to the frequency ujq of the resonator field determines the na- 
ture of the dressed medium. The critical condition of Eq. < fT~8b 
indicates the competition between the population pumping to 
the excited level and the spontaneous relaxations to the ground 
state. The transparency effect is present only when the coher- 
ent pumping is sufficiently strong to overcome the relaxation; 
otherwise, the probe signal is trapped and the dressed medium 
becomes absorbing. 

While the imaginary part of the susceptibility is an even 
function of A, the real part is odd and one-order higher in A 
than the imaginary part. The dispersion spectrum similarly ad- 
mits multiple local extrema, though the dispersion is always 
zero at A = for any value of qubit level spacing uj q be- 
cause of the odd symmetry about the origin. However, when 
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sweeping us q across the resonance point cjq/2, the spectrum is 
inverted, i.e. the dispersion is switched from positive to nega- 
tive or vice versa. 



B. Near-degeneracy and the switching from closed to open 
transitions 

The phenomena discussed in the last subsection are some- 
what analogous to those presented in Refs. lfl3l [3, so 
here we give a physical interpretation of the tuning from 
transparency to absorption using the terminology of atomic 
physics. 

The three dressed states |^t )i |^o) an d l/^i) that we have se- 
lected as the basis of the three-level system have tunable level 
spacings based on the transformation angles 6q and Q\ which 
are defined by the detuning between the parameters ui q and 
wo- After coupling to the strong control field, the two dressed 
states |^o) and |yUi) will have line-broadening. The metastable 
state |fo) in particular might be so broadened that it overlaps 
with the ground state \/j,q). Such a case most easily occurs 
at resonance with u) q = ujq/2, where the spacing between the 
lower two states \fio) and |i/o) is minimized to 2ij (Cf. Eq. ©). 

When the overlap occurs, the three-level system becomes 
effectively two-level. The lower two levels |/io) and \i/q} be- 
come quasi- or near-degenerate; more precisely, they degener- 
ate into a single ground state and differ from each other only 
as hyperfine levels of the common ground state. The minimal 
amplitude il c of the coupling needed for overlapping can be 
roughly estimated using a first-order perturbative expansion. 
Considering the spacing between the metastable state |^o) and 
the excited state to be E^ — Eq = uj — (\/2 + 1)77, we 
find the level shift of j^o) up to first order to be 



u 



{V2 + 1W 



After equating the above to 2rj, the amplitude needed can be 
obtained as 



a. 



2rj 



UJQ - (V2+ 1)?7 



(20) 



Comparing Eq. (T20b with Eq. (fT8l l. we can observe that, 
in order to exhibit electromagnetically induced transparency, 
the coupling amplitude Q c would have to be within a range 
such that it can simultaneously overcome the spontaneous re- 
laxation and yet prevent the system from being degenerate, 
in addition to the requirement that Lo q be outside the critical 
range indicated by Eq. dT9l . Without this range, the dressed 
medium becomes unresponsive to the incident probe signal 
and the absorption spectrum becomes flat. 

If we draw the analogy of the lower dressed levels to the 
Zeeman sublevels of the degenerate ground state lfl3ll . then 
|/Ko) and \vq) can be deemed sharing the same non-zero "an- 
gular momentum." The remaining factor for deciding the 
dressed medium to be transparent or absorbing depends on 
whether the transition between the two sublevels is open or 
closed; in other words, whether the transitions within the 



three-level system are non-cyclic (A-type) or cyclic (A-type). 
Different from the usual case, where this factor is determined 
by the hyperfine structure of a particular atom, the dressed 
qubit system we discuss here has this factor effectively deter- 
mined by the dressed relaxation rates 7 M i and 7 Mi , . 

From Eqs. fl!6dt -( [T6fl ), we observe that unlike the usual 
multi-level SQUID systems, the relaxation rates of the dressed 
qubit are tunable through the transformation angles 9q and 6±. 
In addition, the relaxation rate 7 M „ between the lower levels is 
actually derived from the dephasing rate r v of the undressed 
qubit. The type of transitions in the dressed qubit thus be- 
comes exploitable since the dephasing rates of the supercon- 
ducting qubits are in general much greater than their relax- 
ation rates n and the different dependences on the transfor- 
mation angles control whether > |7 M i| or < 

When the qubit is closely resonant with the CPW resonator, 
or precisely, oj q is within the critical range between Ac,- and 
Ac,+, the magnitude of 7 M „ is greater, under which the flip- 
ping processes of the populations between the ground and the 
metastable states dominate over the excitation process from 
the ground state to the excited state. The transition between 
the ground and the metastable states is thus closed, effectively 
degenerating the two levels and making the three-level system 
operate in a A-type setting. On the other hand, when the qubit 
is off-resonant with the CPW resonator, the opposite condi- 
tion 7^1 > 7^ is met and the lower two levels become suf- 
ficiently non-degenerate that the transitions among the three 
levels cannot be considered cyclic. The system then operates 
in the usual A-type setting that electromagnetically induced 
transparency can take place. 



VII. NUMERICAL ANALYSIS 

We now study the two parts of the susceptibility by consid- 
ering experimentally accessible parameters. We first examine 
the dressed charge qubit, which was experimentally realized 
in Refs. I25ll33ll . We give the variation of the susceptibility 
against multiple parameters and identity the effective ranges 
of the qubit level spacing for the EIT and EIA. The cases for 
phase and flux qubits are discussed later to show that the ar- 
guments for charge qubits can be applied to other qubits, for 
tuning the susceptibility to different operating regimes. 



A. Charge qubit 

Without loss of generality, we now assume that the res- 
onator frequency is fixed at, e.g., ljJq/2tt = 7 GHz. We 
consider the charge qubit model with the following param- 
eters: the qubit has a junction energy Ej/2it = 2.6 GHz 
and a charge energy in the GHz range that we use to tune 
the qubit level spacing by varying the gate reduced charge 
number n g ; the coupling coefficient between them is assumed 
to be t]/2tt = 100 MHz. The undressed relaxation and de- 
phasing times of the qubit are taken as 1/n = 0.7 /is and 
l/ r 2 = 48 ns, respectively. The operating temperature is as- 
sumed to be 20 mK. 
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FIG. 4: (Color online) Real part x' (blue) and imaginary part \" 
(red or yellow) of the normalized susceptibility of the dressed charge 
qubit versus the detuning A and the qubit level spacing uj q . The de- 
tuning abscissa is normalized with respect to A' /2tt = 25 MHz; the 
qubit-level-spacing abscissa is normalized with respect to ui' q /2ir = 
3.5 GHz. This plot indicates the dispersion and absorption spectra of 
a probe signal over different operating ranges. 



Figure |4] plots both the real \' and the imaginary \" parts 
of the susceptibility \ m normalized units as a function of the 
normalized probe-signal detuning A/ A' (with respect to the 
normalizing constant A'/27r = 25 MHz) and the normalized 
qubit level spacing io q /uj q (with respect to the normalizing 
constant L>j' q /2ir = 3.5 GHz) near the resonance range of the 
dressed qubit. The strongly coupled control field can achieve 
a coupling amplitude on the order of 10 MHz l25ll ; the plot 
here assumes a value of fl c /2Tr = 12 MHz. 

For the imaginary part \" (warm colored: red or yellow), 
we are able to observe the absorption peaking at the resonant 
frequency uj q = ujq/2 = 2tt x 3.5 GHz and its immediate 
falloff when uj q is tuned slightly off-resonance. Along with 
the attenuation of the magnitude is the symmetric splitting of 
the peak about the zero detuning point A = at either off- 
resonant side. The zero-detuning point itself falls from its 
maximum value to its local minimum along the path of off- 
resonance, indicating the switching of the dressed medium 
from "maximally absorbing" to being transparent to the probe 
signal. 

The switching phenomenon is better illustrated in Fig. [5ja) 
where the imaginary part \" of the susceptibility is plotted 
for various values of the spacing u g . We note that the suscep- 
tibility obtains a maximum value and a minimum half-width 
with a Lorentzian shape when the dressed medium is resonant 
(the thickened curve for u> q /2-K = 3.5 GHz). Following the 
detuning between the qubit and the CPW resonator, the half- 
width starts to spread out while the peak starts to dent. The 
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FIG. 5: (Color online) Normalized susceptibility spectra at typical 
values of the qubit level spacing ui q plotted versus the scaled detun- 
ing A/ A': (a) the absorption x" = I m [x] an d (b) the dispersion 
x' = Re[x]. The curves that correspond to the resonant frequency 
cjo/2 = 2-7T x 3.5 GHz are thickened and all the numbers indicate the 
value of ui q taken (in unit of GHz) in both (a) and (b). Note that in (a) 
a dip appears at the center of the curve about uj q /2n ~ 3.47 GHz, 
which indicates the switching of the dressed qubit from being ab- 
sorptive (> 3.47 GHz) to being transparent (< 3.47 GHz). In (b), 
the curves that correspond to u q below are solid and those 

above uio/2 are dashed. Note that the solid and dashed curves with 
equal distance from the resonance wo /2 are symmetric counterparts 
of each other. Their roles for positive and negative dispersion across 
the range of the detuning A are exchanged above and below the res- 
onance frequency. 
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FIG. 6: (Color online) Plot of the normalized absorption spectrum 
x" = Im[x] versus the normalized detuning A/A' and the coupling 
amplitude Q, c to the control field. 

switching occurs at the critical value Xc/2tt = 3.40 GHz, 
according to the approximate analytical solution of Eq. dT9b , 
where the qubit-to-resonator coupling 77 is not very large. A 
numerical estimate gives Xg/2-k 3.47 GHz, as can be seen 
from Fig. |5ja). 

The plot of Fig.|4]also shows that the falloff of the absorp- 
tion is even-symmetric about the resonance point uj q = (Jo/2, 
due to the fact that x" is an even function of the dressed re- 
laxation rates 7^1 and 7 M „, which in turn are odd functions of 
the qubit spacing uj q . This symmetry verifies the exchange- 
able roles of the near-degenerate levels |/^o) and \vq): when 
u! q < wo/2, I Mo) nas the lower eigenenergy of the dressed 
states and is regarded as the ground state that couples to the 
probe signal, while \uq) is regarded as the metastable state; 
when ujq > luq/2, \fio) becomes the metastable level and \uq) 
becomes the ground state. 

In Fig. [4] for the real part \' (blue colored), we can see the 
dispersion spectrum having an odd symmetry around the zero 
probe detuning point A = 0, i.e. x'( w <?)|a=o = 0, through- 
out the range of the qubit spacing uj q . Unlike x"< this real 
part x' of the susceptibility is an odd function of the dressed 
relaxation rates, resulting in its odd symmetry of u q around 
the resonance point ui q = ujq/2, which is better illustrated 
in Fig. Ob) . There, the dispersion spectrum for various val- 
ues of w q is plotted versus the detuning A. The one curve at 
resonance (u! q /2n = 3.5 GHz) is thickened and the ones that 
have their symmetric counterparts below resonance are shown 
as dashed curves, from which we can observe the switching 
from positive dispersion to negative dispersion of the dressed 
medium discussed in Sec. I VII 

In addition to the spacing ui q , tuning the dressed qubit from 
transparency to absorption also depends on the dressed sys- 
tem's coupling strength fl c to the control field since this cou- 
pling strength directly affects the absorption (Cf. Eq. (fT5]l). 
Its influence is illustrated in Fig. [6] where the absorption is 



plotted against A/ A' (over the same range as in Fig. [4]i and 
O c (from 4.5 MHz to 45 MHz) while the qubit spacing is held 
at a typical value u q /2Tr = 3.4 GHz. We can notice the sin- 
gle peaking at the lower end of the coupling to the control 
field, where the dressed medium is weakly driven by the con- 
trol field and exhibits a population trapping of the probe sig- 
nal. Towards higher values of the coupling, the excited state 
of the dressed three-level system is sufficiently detuned from 
the probe signal that it starts to exhibit the transparency ef- 
fect. Then similar to ordinary A-type atoms, the dressed qubit 
has its twin absorption peaks further apart when the coupling 
strength is increased. 



B. Phase qubits 

We now examine our general theory of the tunable trans- 
parency and absorption effects on other superconducting 
quantum circuit systems. 

For a phase qubit, we adopt the experimental parameters 
of Ref. DOIl : the CPW resonator has a frequency u}o/2n = 
6.57 GHz; the coupling strength between the resonator and 
the qubit is fixed at i]/2ir = 19 MHz; the undressed relaxation 
and dephasing times of the qubit are, respectively, 650 ns and 
150 ns. The coupling strength to the control signal is assumed 
to be Q, c /2n = 3.85 MHz. The operation temperature is held 
at 25 mK. The normalized absorption spectrum is plotted as a 
contour plot versus the normalized qubit level spacing uj q /uj q 
with respect to the charge qubit case (uj q from 3.2 GHz to 
3.37 GHz; normalizing constant uj q the same as that of the 
charge qubit) and the normalized probe detuning A/ A' (A 
from -6.5 MHz to 6.5 MHz; normalizing constant A' the same 
as that of the charge qubit) in Fig. [3a). 

Similar to the case of a charge qubit in the last subsec- 
tion, the peak at the center of the contour falls off symmet- 
rically as the qubit spacing is tuned off the resonant frequency 
ui q /2ir = 3.285 GHz. The peak also splits symmetrically 
around the zero detuning A = 0, and enters from the ab- 
sorption region into the transparency region at the critical fre- 
quency Xc/2n = 3.266 GHz according to Eq. (fl9l . The half- 
widths of the absorption peaks decrease along with the fall off 
of the magnitude. 



C. Flux qubits 

For a flux qubit, we adopt the experimental parameters 
of Ref. J3lll : the CPW resonator has oscillating frequency 
u!o/2ir — 9.907 GHz; the qubit-resonator coupling strength 
is T]/2ir 100 MHz; the undressed relaxation and dephas- 
ing times of the qubit are, respectively, 1.9 /is and 1 /js 
113911 . The coupling strength to the control signal is set to 
n c /27r = 0.63 MHz. The operation temperature is 50 mK. 
The normalized absorption spectrum of the dressed flux qubit 
is also plotted as a contour plot versus the normalized qubit 
level spacing u) q /uj q with respect to the charge qubit case (w g 
from 4.7 GHz to 5.2 GHz; normalizing constant tu' q the same 
as that of the charge qubit) and the normalized probe detuning 



12 



(a) 



curs at the critical frequency Ac/27r = 4.854 GHz. In this 
case, the absorption peaks also split out slower and have nar- 
rower half- widths. 
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FIG. 7: (Color online) Contour plots of the normalized absorption 
spectra \" versus the normalized detuning A/ A' on the vertical axis 
and the normalized qubit level spacing u] q /uj' q on the horizontal axis: 
(a) for a phase qubit and (b) for a flux qubit. The red end of the color 
spectrum (i.e., the centers of both contour plots) indicates peak val- 
ues of the absorption, where the qubit is maximally dressed with a 
rotation angle 9q = tt/2. The deep blue end of the color spectrum 
for A/ A' = (i.e., the middle sections of the left and right bor- 
ders of both plots) indicates the minimum values of the absorption, 
where the qubit is dressed with a rotation angle 8o < tt/2. Namely, 
maximum transmission due to EIT occurs at those blue ends. 



A/ A' (A from -1.5 MHz to 1.5 MHz; normalizing constant 
A' the same as that of the charge qubit) in Fig.|7tb). 

Similar operating regions of transparency and absorption 
can be observed in this flux qubit system compared to the 
other qubit systems discussed before, except for the scale of 
variations. For instance, comparing Figs. [7ja) and (b), we see 
that the fall-off of the magnitude along the zero detuning is 
relatively slower in the flux qubit case due to the slower de- 
phasing time of the undressed flux qubit. The switching oc- 



VIII. CONCLUSION 

We have proposed a method to realize the effects of both 
EIT and EIA on superconducting quantum circuits using 
dressed states derived from the coupling between an arbitrary 
type of two-level superconducting junction qubits and a copla- 
nar waveguide resonator. The use of dressed states alleviates 
the need to maintain a multi-level structure of the Josephson 
devices and gives rise to tunable relaxation rates between the 
energy levels. The tunable structure of the levels leads to the 
switching between EIT and EIA, which depends on the vari- 
able qubit level spacing and is associated with the open or 
closed transition structure and the hyperfine degeneracy of the 
dressed three-level system. 

Our investigation demonstrates another example of non- 
linear optical phenomena implementable on superconducting 
quantum circuits. We can also see that the special characteris- 
tics of Josephson junction devices, the externally controllable 
Josephson coupling energy in this case, could bring new per- 
spectives to the study of quantum optics where, for example, 
the many parameters are usually fixed for the particular type 
of atom studied and the cavity QED system that surrounds it. 
The switching between EIT and EIA might have important ap- 
plications for the control of superconducting circuits and for 
quantum information transfer in these systems. 
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Appendix A: Condition of local extrema for the susceptibility 

The two roots of Eq. ( fTTI ) that correspond to local extrema 
of the dispersion spectrum must satisfy 



which is equivalent to 

When squaring the two sides, the above inequality implies 

7uiCc + 2 7mi7m"Cc ~ 7ui7L > • 
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For Q > —Jfj.ijfj.iy, with the relaxation rates taking values so the real roots exist under this condition. 
7„i < and > from Eqs. ( |16dt -( |16fT l, we have 

7„i{(7m1+ 2 7^)Cc -7)1} 
> {-(7 M i + Z-ffjuhfjiJuv ~ llu) 
= -7^17^(7/^1 + 7^) 2 > °! 
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